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Abstract 

This paper proposes a gravitodynamic theory because there are similarities be- 
tween gravitational theory and electrodynamics. Based on Einstein's principle of 
equivalence, two coordinate conditions are proposed into the four-dimensional line 
element and transformations. As a consequence, the equation of motion for gravi- 
tational force or inertial force has a form similar to the equation of Lorentz force 
on a charge in electrodynamics. The inertial forces in a uniformly rotating system 
are calculated, which show that the Coriolis force is produced by a magnetic-type 
gravitational field. We have also calculated the Sagnac effect due to the rotation. 
These experimental facts strongly support our proposed coordinate conditions. In 
addition, the gravitodynamic field equations are briefly discussed. Since only four 
gravitational potentials (3 + 1 split) enter the metric tensor, the gravitodynamic 
field equations in "3+1 split" form would be analogous to Maxwell's equations. 
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1. Introduction 

In a recent review paper, Norton^ has shown that the principle of general 
covariance has been disputed as physically vacuous. The debate over this principle 
persists today. In general relativity, space-time coordinates can be arbitrarily 
chosen, I 2 ' 3 ' and gravitational potentials are represented by the metric tensor, g^ v , 
in the four- dimensional (4-D) line element: 

ds 2 = g iiv dx il dx v (1) 

where g^ v has ten independent components, and Greek letters ji and v run from 
to 3. In Eq.(l), ds 2 is an invariant under space-time coordinate transformation. 
However, in order to get unique solutions of the Einstein field equation, certain 
coordinate conditions must be imposed. Ref. [4] is a recent example that discusses 
the choice of coordinate conditions. In this paper, we emphasize that constraints 
on coordinate conditions are imposed by the principle of equivalence. 

Einstein illustrated the principle of equivalence using the following simple ex- 
ample. Let E be an inertial reference system; let another system K be uniformly 
accelerated with respect to E. Then relative to K, all free bodies have equal and 
parallel accelerations. They behave just as if a gravitational field were present and 
K were unaccelerated.t 2 '^ In another words, inertial force is a type of experimen- 
tally confirmed "gravitational force" . 

In order to formulate the above simple equivalent gravitational field suggested 
by Einstein, let us start with an inertial system E. Let X 1, = (cT,X,Y,Z) = (cT, 
R) be the pseudo-Cartesian coordinate in the system E. Then the 4-D line element 
takes the form 

ds 2 = (cdT) 2 - dX 2 - dY 2 - dZ 2 (2) 

We now consider a system K with coordinates a/ 4 = (ct,x,y,z) = (of, r), which 
has a constant acceleration, a, with respect to E. Since the Lorentz transformations 
are not valid between E and K, the following quasi- Galilean transformations are 
usually introduced: ^ 

R = r + iat 2 ; T — t (3) 

A simple calculation for Eqs.(2) and (3) gives 

ds 2 = (1 - (at) 2 /c 2 )(dx°) 2 - 2(af/c) • drdx° - dr 2 (4) 

Comparing (4) with (1), we obtain g OQ = (1 — a 2 t 2 /c 2 ); g Q i = — Oji/c; gij = — 
Latin letters i and j run from 1 to 3. 

We substitute these components of g^ v into the geodesic equation giving 

d 2 x^ dxPdx x 
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In the case of low velocity approximation, Eq.(5) reduces to d 2 r/dt 2 = —a. This 
indicates that a particle with a rest mass m in a uniformly accelerated system 
should experience a uniformly inertial force (— m a). Although this result is well 
known, the above calculation shows that certain constraints on space-time coordi- 
nates and transformations should be imposed by the principle of equivalence. As is 
well known, the inertial force is a type of measurable force. Just like other forces, 
inertial force is a vector, which can be reduced or balanced by other forces. How- 
ever, if arbitrary curved coordinates are chosen, in general, the force components 
cannot be singled out from the expression of Eq.(5). Then the physical meaning 
would be lost in the mathematical maze. This is a major reason whythe debate on 
the principle of general covariance still persists. M Only when the equation of mo- 
tion in a vector form is derived, curvilinear spatial coordinates can be introduced 
for different applications. 

2. Two Proposed Coordinate Conditions 

Based on the principle of equivalence, we propose two coordinate conditions 
into the 4-D line element and transformations in our formulation: 

(i) In any reference system, the 4-D line element takes the standard form 

ds 2 = V {dx°) 2 - 2(* • dr)dx° - dr 2 (6) 

Comparing Eq.(6) with ds 2 = g^dx^dz" , we have g 0O = \I/ , g oi = — <?y = —Sij. 

(ii) Between two reference systems, K and K', the value of ds 2 is an invariant, 
but space-time coordinates take infinitesimal Galilean transformations 

dr = dr - vdt; t' = t (7) 

where v = v(r,t) is the 3-D relative velocity between two systems K and K' . As 
a special case, one reference system can be an inertial frame which has Minkowski 
metrics. 

The physical implication of the above proposed conditions is described by 
Weinberg.® He has emphasized that a local gravitational field has two kinds of 
sources. One is nearby mass; another is all the mass in the universe. It can 
be anticipated that the inertial systems are determined by the mean background 
gravitational field produced by all the matter of the universe. This background 
gravitational field (or simply called ether) provides a base for us to propose the 
above two specific coordinate conditions. 

Landau et al' 7 ' and M0ller® have rewritten the 4-D line element, Eq.(l), in 
another form: 

ds 2 = ( lo dx° - -udx 1 ) 2 - da 2 (8) 
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Comparing Eq.(8) with Eq.(6), 7 D = ^fgZ = ^/W ] ^ = -g oi / JgZ = ^/W Q ; da 2 
= 7 ifc cfafefc^, where 7^ = -# ife + 7^ = 5 ik + *i* fe /* . da 2 is the 3-D measured 
spatial interval in system K. The proper time is defined as dr = ds/c, which is an 
invariant and indicates the measured time interval in system K. For slow motion 
and weak field approximation, it is easily seen that the coordinate time interval 
dfc^dr , and the coordinate spatial interval dr 2 ~da 2 . 

Notice that fy and \I/ are components of metric tensor g^ u in Eq.(6). Although 
they have properties of gravitational potentials, they do not form a 4-D vector. 
From Eq.(6), it is easy to derive the transformations of these potentials 

2 

tt ' = * - \ - * • -; = * + - (9) 

<r c c 

These transformations form a group with the parameters v. 

3. Gravitational Equation of Motion 

In this section, we will rewrite the general equation of motion, Eq.(5), into an 
explicit form in terms of the proposed two conditions. We consider the variational 
principle ^ 

8 J {-m c)ds = 8 j L{x\u\t)dt (10) 
where L^^J) is the Lagrangian. Using Eq. (6), we obtain 



L(x\u\t) = (-m c 2 )^ -2V -u/c-u 2 /c 2 (11) 

where u = dr/dt is the velocity of a particle in system K. 

Let y&o = (1 + 20). Substituting the Lagrangian (11) into the Lagrange 
equation, ' 7 1 

d ,dL, dL 

« ( ^ )= aF (12 > 

we obtain the equation of motion: 

d(mu) T dm . u , , . . 

1 + cV— = m(g + -xh) (13) 
dt dt c 

where g and h are defined as 

g = c\-V<p-- c ^y, h = c 2 ( V x*) (14) 

Here, g is gravitational intensity and h is magnetic-type gravitational intensity. 
Both of them have the units of cm/sec 2 . Therefore, the gravitational force is 
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analogous to the Lorentz force on a charged particle in an electromagnetic field. 
As is discussed in section 1, force is a vector, which can be reduced or balanced by 
other forces. We now have clearly singled the force term out in the right hand side 
of .Eg. (13). However, there are two different terms between the Lorentz formula 
and Eq. (13): (i) a moving mass m = m T, where T = (\l/ — 2\E' • u/c — u 2 / c 2 )~ l l 2 , 
while the electric charge is a constant; (ii) a small additional term on the left hand 
side of Eq. (13), which is related to the changing rate of a moving mass. We stress 
that Eq. (13) is the explicit formulation of the gravitational equation of motion, 
which is also valid to describe the magnetic-type gravitational force as well as 
inertial force in non-inertial systems. 

We now consider the total energy of a moving particle with a velocity u in 
a gravitational field. According to Lagrange's formulation, the total energy of a 
particle is defined as' 7 ' 

E t = ^--u-L (15) 
Substituting the Lagrangian expression (11) into (15), one obtains 

E t = m o rc 2 [(l + 20) -u/c] (16) 
For a weak field, (i,f < 1, an approximation of Eq. (16) is 

Et = m ; c2(1 + 0) (it) 

'l - u 2 /c 2 



Eq. (17) is essentially the same as the formula derived by Landau and LifshitzJ 7 ' 



4. Inertial Force in Accelerated Systems 

In terms of the formulation described above, inertial force is a type of exper- 
imentally confirmed gravitational force. We now calculate two typical cases of 
inertial force. 

(a) . A Uniformly Linear Accelerated System 

In the introduction of this paper, we discussed the inertial force in a uniformly 
linear accelerated system. From Eqs.(4) and (6), we have <fi = — a 2 t 2 /2c 2 , = 
sd/c. By using Eq. (14), we obtain 

g = -c— = -a; h = (18) 

Notice that this equivalent gravitational intensity g in Eq. (18) is produced by the 
3-D gravitational vector potential. 

(b) . A Uniformly Rotating System 
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Suppose a uniformly rotating system has a constant angular velocity O = f2k 
with respect to an inertial system E. Let the infinitesimal Galilean transformation 
be 

dr = dR —(fix R) dT; t— T (19) 

Substituting Eq. (19) into Eq.(2), we obtain the 4-D line element in the rotating 
system K as follows: 

ds 2 = (1 - (ft x r) 2 /c 2 )(dx°) 2 - -(ft x r) • drdx° - dr 2 (20) 

From Eq.(20), we have <p = -(ft x r) 2 /(2c 2 ), * = (SI x r)/c. Using £g.(14), we 
obtain 

g = I y (ft x r) 2 = r2 2 rn; h = yx (f]xr) = 2cQ (21) 

where ry = od + yj. -Eg. (21) gives the centrifugal force, F a = mg = mQ 2 r\\, and the 
Coriolis force, F& = mu x h/c = 2muxft in the rotating system. Therefore, the 
Coriolis force is produced by a magnetic-type gravitational field. M0ller's^ derived 
result is only approximate because he did not explicitly choose the coordinate 
conditions. 

5. Sagnac Effect 

In a sagnac interferometer, the light beam coming from a source is split into two 
sub-beams. One sub-beam circulates a loop in a clockwise direction, and another 
one circulates the same loop in a counter-clockwise direction. When the whole 
interferometer with light source and fringe detector is set in rotation with a rate 
of Q rad/sec, a fringe shift A with respect to the fringe position for the stationary 
interferometer is observed, as first reported by Sagnac in 1913' 8,9 l The formula is 
given as 

A = AQAf/c 2 (22) 

We now derive the above formula simply using the 4-D element in Eq. (20). 
As is well known, the propagation of light has zero value of the 4-D line element, 
that is ds = 0. If a light beam travels along a circle with radius R in a uniformly 
rotating system, Eq. (20) can be simplified as follows: dr =Rd09 o , ft x r = HR O , 
then we have 

ds 2 = (c 2 - (VLRf)dt 2 - 2(ttR 2 )d6dt - R 2 (d6) 2 = (23) 

Solving the quadratic equation (23) for Rd0/dt, we obtain the speed of light in 
the frame of the rotating disk: 
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c e = RdO/dt = c + OR (24) 

The sign of minus or plus in Eq. (24) depends on the direction of the light beam. 
When one beam circulates a loop in a clockwise direction, and another one circu- 
lates in the opposite direction, the time difference is 

2irR 2ixR a „ . 2 

At = — — ~ AnR 2 Q c 2 25 

c-QR c + VlR 

Since the shift of fringes equals At multiplying the frequency of the light, 
Eq.(25) leads to A= 404//c 2 , which is the same as Eq. (22). 

As a matter of fact, the Sagnac effect has clearly shown that (i). one-way speed 
of light is direction-dependent in a rotating system, which is given in Eq. (24); (ii). 
the transformation of space and time coordinates between an inertial system and 
a rotating system obeys the infinitesimal Galilean transformations, instead of the 
Lorentz transformations. 



6. Brief Discussion on the gravitodynamic Field Equations 

In this paper, the term "gravitodynamics" is used as a replacement for the 
term "gravitation" since there are similarities between gravitational theory and 
electrodynamics. For instance, in the systems with weak gravitational fields and 
low mass velocities (u <C c), Einstein's field equations can be written in a form 
analogous to Maxwell's equations: I 10 " 14 ] 



V • g = -4:nGp m (26a) 

vxh = iV ( p m u+-^) (266) 

c c ot 

1 dh 

vxg = "c¥ (26c) 

V • h = (26d) 

In Eqs.(26), 4-D space-time is sliced into 3-D space plus 1-D time (3+1 split). Be- 
cause of the different approach employed to approximate Einstein's field equations, 
the deduced Maxwell-type field equations in Ref. [10-14] have small differences with 
Eqs.(26) in high order terms. For instance, the right hand side term in (26c) is a 
high order term, which is often set to zero. Regardless of these small differences, 
Thome' 10 ! pointed out that there are three major differences between Eqs.(26) and 
Maxwell's equations: (i) minus signs in the source term of Eqs.(26a) and (26b) are 
due to the fact that the gravitational field is attractive rather than repulsive; (ii) a 
factor tV=4 on the right hand side of Eq. (26b) presumably is due to the fact that 
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gravitational potential has ten independent components; (iii) a charge density is 
replaced by mass density p m times Newton's gravitational constant. Furthermore, 
Thorne emphasized the practical importance of such "3 + 1 split" for experimental 
physicists. As a matter of fact, preparations are already underway for experiments 
to search magnetic-type gravitational field. I 10 ' 11 ' 15 ] i n principle, the properties of 
gravitational waves can be derived from the field equation (26). However, deriva- 
tion of the velocity of gravitational waves should require using rigorous forms of 
the gravitodynamic field equations. As an example, it is easily seen that the value 
of N in Eq. (26) would affect the velocity of gravitational waves. 

Following the regular procedures to derive Einstein's gravitational field equa- 
tions, one could apply the variation principle: 

5(S mf + S f ) = 8 J (C mf + C f )dVdx° = (27) 

where dV dx°=dxdydzdx° is 4-D volume element. C m f is the part of Lagrangian 
density that depends on the interaction between particles and gravitodynamic field, 
while Cj is the part that depends on the gravitodynamic field itself. Both of them 
could be derived from Eq.(6) and its metric tensor. 

We emphasize that the metric tensor in our formulation has only four non- 
constant components, \l/ and Therefore the derived gravitodynamic field equa- 
tions would have four component equations (3+1 split), instead of ten equations as 
in general relativity. For this reason, the gravitodynamic field equations for a weak 
field would reduce to a form analogous to Maxwell's equations. Such equations are 
approximately shown in Eqs.(26) with a constant N—l in our formulation. 

Generally speaking, the rigorous "3 + 1 split" form of the gravitodynamic 
field equations should be a set of non-linear differential equations in terms of ^/ Q , 

g and h. This approach deserves further investigation. Furthermore, since 
the gravitational equation of motion and gravitodynamic field equations in our 
formulation can be written in 3-D vector forms, one still has freedom to choose 
3-D curvilinear spatial coordinates in a given reference system. At this stage, the 
process of general covariance has its limited application. 
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Appendix 

A Derivation of the Gravitodynamical Field Equations 

In the following, a possible form of the gravitodynamical field equations will be 
derived from Eq. (27). First, C m f can be obtained from Eq. (11). 

Cmf = (p m c 2 ) v /^ D -2* -u/c-u 2 /c 2 (28) 

For the weak field and low mass velocities, Eq. (28) can be approximately rewritten 

as 

Cmf = ~pm(-C 2 + U 2 /2 - C 2 + C (* • U )) (29) 

Eq. (29) is quite similar to the £ m /(EM) in electrodynamics. 

Secondly, regarding the function of Cf, we need to consider the differences 
between the "actual" gravitational field and the "equivalent" gravitational field 
produced by acceleration in a non-inertial system. Here we only intend to derive 
the field equations for the "actual" gravitational field produced by gravitational 
mass. For this purpose, we should study the transformations of field intensities 
here. 

Start with Eq.(9), the transformations of potentials and St. It can be 
rewritten as 

4>' = <l>-~---9- t *' = * + - (30) 
2 c 2 c c 

where v is the relative velocity of system K' with respect to system K. Since 
we only consider the actual gravitational field in both systems, a constant v is 
assumed here. Using the definition of g and h in Eq. (14), Eq. (30) leads to 

g' = g + (3xh; h' = h (31) 

Therefore, the magnetic-type gravitational intensity is an invariant 3-D vector. 
From Eqs. (30) and (31), it is easy to prove another invariant 3-D vector g, which 
is defined as 

g = g- *xh = g'-*'xh' (32) 

In terms of the above invariant, various kinds of Cf can be introduced. Com- 
paring the function of £/(EM) in electrodynamics, we choose the following Cf for 
discussion: 

C f = -(8nG)- l (g-g-h-h) (33) 
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Substituting Eq. (32) into (33), and neglecting the high order terms, we obtain 

C f = -(87rG)" 1 (g • g - 2g • (* x h) - h • h) (34) 

Substituting the expressions of C m f in Eq. (29) and C m f m Eq.(34) into the 
Eq. (27), and following the same procedures to derive the field equations in electro- 
dynamics [7], we can obtain two gravitodynamic field equations: 

V • (g - * x h) = -AnGp m (35a) 

Vx(h-*xg) = -— /) m u + -ig-*xh) (356) 

c cot 

When comparing Eq.(35a, b) with Maxwell's equations, it is easy to see that there 
are additional terms, \P x h and \I/ x g, in Eq.(35a, b). 

Using the mathematical properties of the operator y, other two gravitational 
field equations can be easily derived from Eq.(14): 

VXg = "cI (35C) 

V • h = (35<Z) 

Furthermore, from Eq. (35a) and Eq. (35b), the equation of continuity can be ob- 
tained 

^f + V-(AnU) = (36) 
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